Measurement of ift(2S) decays to baryon pairs 
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A sample of 3.95M ip(2S) decays registered in the BES detector are used to study final states 
containing pairs of octet and decuplet baryons. We report branching fractions for ip(2S) — > pp, 
AA, E°E°, H"H + , A++A", E + (1385)E _ (1385), H°(1530)H°(1530), and n~Ti + . These results 
are compared to expectations based on the SU(3)-&&vot symmetry, factorization, and perturbative 
QCD. 
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I. INTRODUCTION 



In the quarkonium model, the tp(2S) is the first radial excitation of the 3 S cc bound state. As such, its properties 
are expected to be relatively straight-forward to understand, at least in terms of those of the J/tp ground state. 
Somewhat surprisingly, these expectations do not always hold. In particular, there is a rather dramatic anomaly 
associated with the 4'(2S). 

The major puzzle in hadronic tp decays is the large discrepancy between the decay widths for J /ip(lS) — ► pir and 
K*K and the corresponding widths for ip(2S) decays. These modes are expected to proceed via ip — > ggg, with widths 
that are proportional to the square of the cc wave function at the origin, which is well determined from dilepton decays. 
The predicted ratio of branching fractions from factorization is: 



B^(2S) -> Xhad) 

B(J/^ -» X had ) 



a.(il>(2S)) 



a s (J/^j) 
0.116 ±0.022 



B{tp(2S) -> e+e") 
B(J/ip -> e+e-) 



where Xhad designates any exclusive hadronic decay channel. The terms come in from the three gluon widths. ]lj] 
Experimentally, the ip(2S) — > pir and K*K are reduced by over a factor of twenty from these expectations Q. This 
anomaly calls into question the underlying assumption behind the theoretical predictions: that the ip(2S) is a pure 
cc state. 



A. 4>{2S) -> B t B t 

In the context of flavor SU(3), a pure cc state is a flavor singlet and, in the limit of SU(3) flavor symmetry, the 
phase-space-corrected reduced branching fractions to any baryon octet pair, |Mj| , where 

2 = BM2S) 

{p* is the momentum of the baryon in the tp(2S) rest frame), should be the same for every octet baryon, Bi. Deviations 
from this rule could indicate a non-cc component of the charmonium wave function. The reduced branching fractions 
for J/ip — * BiBi decays are shown in Fig. l|. The SU(3) relation works reasonably well, although there may be some 
increase for the pp mode. 




FIG. 1. PDG values for the reduced branching fractions \Mi\ 2 = B(J/ip — > BiBi) /(ivp* / \/s) for J/ip —> BiBi, where 
B t G {p,A,£°,E-,A++,£ ± (1385)}. 
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This relation has not been tested for the ip(2S), where the only relevant mode that has been measured is pp, and 
that with rather poor precision ||, @j. 

There are very few direct calculations of the decay of charmonium to baryonic final states. One of the most 
comprehensive is the perturbative analysis by Bolz and Kroll A comparison to this analysis will be discussed 
later. 



II. THIS EXPERIMENT 

We report results of measurements of the branching fractions for ip(2S) — > Sj-Bj, where Bi £ {p, A, £°, S~, 
A ++ , £+(1385), S°(1530), f2~} using a sample of 3.95 x 10 6 ip(2S) events produced via e + e~ annihilations at the 
BEPC collider and observed by the BEijing Spectrometer (BES). The data represents a total integrated luminosity 
ofwGJpb- 1 . 

The Beijing Electron Spectrometer, BES, is a conventional cylindrical magnetic spectrometer, coaxial with the 
BEPC colliding e + e~ beams ||. A four-layer central drift chamber (CDC) surrounding the beampipe provides trigger 
information. Outside the CDC, a forty-layer main drift chamber (MDC) provides tracking and energy-loss (dE/dx) 
information on charged tracks over 85% of the total solid angle. The momentum resolution is a p /p = 0.017^1 +P 2 (p 
in GeV/c), and the dE/dx resolution for hadron tracks is ~ 11%. An array of 48 scintillation counters surrounding the 
MDC provides time-of-flight (TOF) information of charged tracks with a resolution of « 450 ps for hadrons. Outside 
the TOF system, a 12 radiation length, lead-gas barrel shower counter (BSC), operating in self-quenching streamer 
mode, measures the energies of electrons and photons over ps 80% of the total solid angle. The energy resolution 
is cte/E = 0.22/v^E (E in GeV), and the spatial resolutions are = 4.5 mrad and a z = 4 cm. Surrounding the 
BSC is a solenoidal magnet that provides a 0.4 Tesla magnetic field in the central tracking region of the detector. 
Three double layers of planar counters instrument the magnet flux return (MUID) and are used to identify muons of 
momentum greater than 0.5 GeV/c. Endcap time-of-flight and shower counters extend coverage to the forward and 
backward regions. 

III. BARYON OCTET 

A. ip(2S) -> pp 

The experimental signature for the decay ip(2S) — > pp is two back-to-back, oppositely charged tracks each with a 
momentum of 1.586 GeV/c. The proton typically deposits one-half or less of its 0.91 GeV kinetic energy in the BSC; 
the antiproton undergoes an annihilation process in the BSC approximately half the time, producing a large shower. 

Major potential backgrounds are: ip(2S) — > K + K ~ , 7r + 7r~ , fi + fJ,~, and e + e~ . Each of these modes has a momentum 
at least 190 MeV/c greater than that of the pp channel. 

We select events with two and only two well reconstructed, oppositely charged tracks with good time of flight 
information, and which are not identified as muons by the muon system. Also |cos0| must be less than 0.6 for both 
tracks to ensure that they occur within the fiducial volume covered by the muon system. Candidate pp pairs are 
required to be within 1.8 degrees of collinear. 

The shower counter energy deposition as a function of momentum for positively charged tracks is shown in Figure |^. 
The faint cluster near p — 1.6 GeV/c, Esc = 0.3 GeV is the proton signal. The other features on the graph are 
due to Bhabhas (large concentration at p w 1.8 GeV/c, Esc ~ 1-5 GeV), muons (vertical stripe at p « 1.8 GeV/c, 
E S c < 1 GeV) and radiative Bhabhas (trailing cluster at p < 1.6 GeV/c, E S c ~ p)- To remove these backgrounds, 
a cut is made at Esc < 0.7 GeV/c. In addition, the shower counter has a number of support ribs which are dead 
regions, thus degrading the energy measurement. Tracks which enter these regions are removed from consideration. 
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FIG. 2. Shower counter energy vs. momentum for positively charged tracks. Signal is expected near a momentum of 
1.6 GeV/c, and energy of 300 MeV. 



An additional handle on the identification of protons is gained from the dE/dx system. Figure |^ shows the dE/dx 
particle ID results for candidate events that pass the above cuts. Units are \ — \dE / dx mcas — dE/dx cxp \ /a, where a 
is the resolution of the particle ID system. The vertical axis is for the p hypothesis, and the horizontal refers to the p 
hypothesis. The cluster near (0,0) contains true pp events, and the cluster near (5,5) is a mixture of event types such 
as radiative Bhabhas and ip' — > ee. A cut is made on the combined \i (Xp + xi) 5 < 3. 
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FIG. 3. Distribution of \v f° r antiproton candidates versus Xy f° r proton candidates, where \v — \dE / '<Ar moas — dE/dx exp \ 
and is calculated assuming the track to be a proton (or anti proton). The signal is the cluster near (0,0). The cut made is 

(Xp + Xf)* <3. 

The weighted average momentum spectrum of the remaining candidate events is shown in Figure ^. By weighted 
average we mean that the track parameters of the positive and negative tracks (curvature and dip-angle) are averaged 
together and then combined to form a momentum. This spectrum in Figure Q is fit to a gaussian plus a quadratic 
background function, with the centroid of the gaussian fixed to the theoretic momentum of the protons, 1.586 GeV/c. 
The width and height are allowed to vary. From the fit, N p p = 201 ± 14 ± 20. Here and below, the first error is 
statistical and the second is systematic, in this case the error on the fit. 
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FIG. 4. Weighted average momentum of pp pairs, fit to a gaussian plus a quadratic. 



B. iP(2S) -> AA 

The decays ip(2S) — > AA produce two back-to-back As, each with momentum 1.467 GeV/c. We only consider 
events where both As decay to the charged pir final state. The final states of interest are thus, i/j(2S) — > ir + ir~pp, 
where the p7r~ and pn + originate from well separated decay vertices. The decay kinematics are such that the proton 
(antiproton) is always the highest momentum positive (negative) track in the event. 

We select events with four and only four well reconstructed tracks with a zero net charge, and in the fiducial region 
covered by the drift chamber, |cos0| < 0.80. Events which pass these cuts are processed through a detached vertex 
finding algorithm, and subjected to a 5-C kinematic fit to ppn^Tt" , with M p7T - = Mp^+. The 84 events which pass 
this fit with a confidence level of more than 1%, and have M p7T < 1.15 GeV/c 2 are shown in Figure [j| Extrapolating 
the two events in the region above 1.13 GeV/c 2 and below 1.15 GeV/c 2 to the area under the mass peak, we find that 
there are four background events in the plot. We conservatively assign this number a 100% error and determine iVVr 
tobe80±9±4. 
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FIG. 5. Distribution of pir invariant masses from a kinematic fit to ip(2S) — > pp-K + ir~ , M p7T - = M- n+ . Upper figure is full 
range of M - , lower figure is expanded near the signal peak at 1.11 GeV/c 2 . 



C. i>(2S) -> S°E° 

The E° hyperons from ip(2S) — > decay promptly via S° — > 7A. We consider only those decays where the 

daughter As decay via the charged pir mode. The experimental signature is thus, ip(^S) — > pj)7r + 7r~77, where the 
p7r _ and ]wr + originate from A hyperons with well separated decay vertices. In addition, there are two photons in 
the energy range 27 < E 1 < 202 MeV. As in the case for tp(2S) — > AA, the proton (antiproton) is always the highest 
momentum positive (negative) track in the final state. 

We extract ip(2S) -> event candidates (E° -» A7, A — > pir ) using the same selection criteria as used for the 

AA mode with the additional requirement that there be two or more isolated clusters in the BSC with energy greater 
than 60 MeV, and within region |cos6*| < 0.75. By "isolated" we mean more than 12.8° (cos#i so i < .975) away from 
each of the charged tracks. 

Both p7r pairs in the surviving events are processed through a displaced vertex-finding algorithm and the event 
is then subjected to a five-constraint kinematic to the hypothesis ip(2S) — > p|wr + 7r - 77, with the beam constraint 
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M pw - 7 = M™.+ 7 . Here the highest momentum positive (negative) track is classified as the proton (antiproton) . For 
events with more than two 7 candidates, the fit is applied for each possible combination. 

Events which pass the kinematic fit with a confidence level greater than 1%, and M p7ri < 1.3 GeV/c 2 are shown in 
Figure |[ We fit this spectrum to a single gaussian plus a linear background with the peak fixed to the mass of the 
S°, 1.192 GeV/c 2 . From the fit, N^ = 8 ± 3 ± 2. 
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FIG. 6. Distribution of p-iry invariant masses from a kinematic fit to ip(2S) — ► ppix + Tx 77, M 



Events with 



pnj masses below 1.3 GeV/c 2 are fit to a gaussian signal plus a linear background. There are 8 ± 3 ± 2 events in the E° peak. 



D. if)(2S) -> 



The S _ hyperon from ip(2S) — > decays via S _ — > 7r~A. We consider only those decays where the daughter 

As decay via the charged pir mode. The experimental signature is thus ip(2S) — > p|57r + 7r~7r + 7r~ where one each of 
the £»7r~ and p7r + combinations originate from A hyperons with well separated decay vertices. As in the case for 
ip(2S) — > AA and the proton (antiproton) is always the highest momentum positive (negative) track in the 

final state. 

We select events with six and only six well reconstructed tracks with zero net charge, and in the fiducial region 
covered by the drift chamber, |cos 6*| < 0.80. Each of the four possible pir~ and pir + combinations are sent through 
a displaced vertex-finding algorithm and subsequently subjected to a five-constraint kinematic fit to the hypothesis 
ip(2S) — > pp-K + 7r~ 7r + ir~ , with the beam constraint M p7T -^- = Mp n + n +. 

Events which pass the fit with a confidence level greater than 1% are examined further. We additionally require 
that the pir combinations have a mass within 10 MeV/c 2 of the A and that the mass of the AA candidate is more than 
20 MeV/c 2 away from the J/ij; in order to reduce background from the cascade decay ip(2S) — * J/i/jtttv, J/ip — * AA. 

The M p7r -n- spectrum of events which remain after the above cuts is plotted in Figure 0. There are 12 ± 3.4 
events in the S~ peak. Averaging the five events outside the peak region over the entire plot and multiplying by the 
width of the signal gives 0.15 background events. A conservative error of 100 percent is applied, giving 12 ± 3.4 ± 0.2 
ijj(2S) — > S~S + events detected. 
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FIG. 7. Distribution of piv~iT~ masses from a kinematic fit to ip{2S) — > ppir + tt~ tt + tt~ , M p1T - n - = M- n + n +. There are 
12 ± 3.4 ± 0.2 events in the H~ peak. 

IV. BARYON DECUPLET 

A. ip(2S) -> A++A~ 

The decay ^{2S) ->■ A++A — produces back-to-back A++ and A — . As the A++ is a broad (111 MeV/c 2 ) 
resonance, the primary hyperons do not have well-defined momenta, in contrast to the octet cases above. We select 
events where both A ++ and A~~ decay to pit [i3(A ++ — > pit) > 99%]. The final state is ip(2S) — > pp7r + 7r~. 

We select events with four and only four well reconstructed tracks with a zero net charge, and in the fiducial region 
covered by the drift chamber, |cos#| < 0.80. The surviving events are processed through a four-constraint kinematic 
fit to the hypothesis ip(2S) — * ppir + ir~. Events which pass with a confidence level greater than 1% are examined 
further. 

Figure |8| shows the invariant mass distribution of the pp pair in events which pass the fit. There is a clear peak 
in the J/ip mass region coming from the cascade decay ifi(2S) — > J/ipir + Tr~ , J/iJj — > pp; we remove this by making a 
60 MeV/c 2 cut around the J/ip. Fi gure H shows the invariant mass distribution for pn containing a peak at the A 
mass. We remove the AA background by requiring the pir~ and pn + masses to be greater than 1.15 GeV/c 2 . 
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FIG. 8. Distribution of n + ir recoil masses in the A ++ analysis. A cut is made at |M P p — JWj/^l > 60 MeV/c 2 to remove 
J/tp contamination. 
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FIG. 9. Distribution of piv masses in the A ++ analysis. A cut is made at M - > 1.14 GeV/c 2 to remove AA background. 



Events which pass all above cuts are fit to a spin-1 Breit-Wigner plus a 4-body phase-s pac e background histogram. 
The width and centroid of the signal spectrum are fixed to the PDG M values. Figure]!^ shows the output of the 
fit; there are 849 total events in the plot. The fit parameter varied is the relative proportions of the phase space 



background and the Breit-Wigner signal to the total number of events in the plot. N 



- is 157 ±13 ±34. 
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FIG. 10. Distribution of pit masses from a kinematic fit to ip(2S) — > ppn n~ . Data is fit to a spin-1 Breit-Wigner 
plus a 4-body phase-space background (from Monte Carlo). Black boxes with error bars are data, smooth curve is the 
spin-1 Breit-Wigner fit result, and histogram is the final fit to background plus Breit-Wigner, binned to match the data. 
n a++a— = 157 ±13 ±34. 



B. iP(2S) -> E + (1385)S-(1385) 

The hyperons from ip(2S) — > E* + E*~ decay via S + (1385) — > Att + 88% of the time. We consider only those decays 
where the daughter As decay via the charged pir mode. The experimental signature is ij){2S) — > pp2{ir +, K~) where 
one each of the p-n~ (j>tt + ) candidates is consistent with being from the decay of a A (A). 

We select events with six and only six well reconstructed tracks with a zero net charge, and in the fiducial region 
covered by the drift chamber, |cos 9\ < 0.80. We kinematically fit the 36 possible charge combinations of (H — ) + (— +)— , 
running the (H — )/( — h) candidates through a displaced vertex finding algorithm, to p7r~7r + p7r + 7r~. No constraints 
are placed on the (+—)/( — b) candidates. 

Figure |ll| shows that the fit mass of the daughter As from the decay of the primary S* + is well defined and centered 
at the A mass. We make a loose cut of 15 MeV/c 2 on the A and A resonances, indicated by the arrows on the plot. 
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FIG. 11. Distribution of M pw - from a fit to V(2S) -> (p7r~)7r + (p7r + )7r _ , showing the A peak in E(1385) + E(1385) _ candidate 
events. 

As shown in Figure the mass recoiling against the orphan 7r + 7r~ pair is dominated by a peak at the J /ip mass, 
indicating contamination of i/>(25 f ) — * tt + tt~ J/ip, J/ip — ► AA. We therefore remove events with a n + n~ recoil mass 
within 30 MeV/c 2 of the J/ip. 
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FIG. 12. Distribution of tt + tt recoil masses from a kinematic fit to ip(2S) — > AA-k + tt , showing large J /ip contamination 
in E(1385) + E(1385)" candidate events. 

To determine iV s ,. + ^,_, we constrain the J)7r + 7r~ combination candidate) to be within 107.4 MeV/c 2 (3 x Y) 

of the nominal PDG value in order to enhance the S* + signal. Events which pass the above cuts are fit to a 
Breit-Wigner with a constant background, with the mass and width fixed to the PDG values (M = 1382.8 GeV/c 2 , 
r = 35.8 MeV/c 2 ). This fit is shown in Figure [If; from the fit, N^ t+ ^,_ = 13.8 ± 3.7 ± 2.7. 
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FIG. 13. Distribution of M pn - n - from a kinematic fit to ip(2S) —> 
Breit-Wigner constrained to the nominal E(1385) + mass and width. 



(pn )n (pir + )ir + . The histogram is an unbinned fit to a 



C. V(25) -> H°(1530)H°(1530) 



In the decay i/j^S) ~> the S*s are produced back to back in the i^(2S) rest frame. The do mina nt decay 

mode of S* baryons is S* — > S~7r+, with a branching fraction of 0.66. |?J The S~ decays as in Section HID to An~ , 
and the A decays to pir~ . 

We select events with eight and only eight well reconstructed tracks with zero net charge, and in the fiducial region 
covered by the drift chamber, |cos(9| < 0.80. Remaining events are subjected to a 4-constraint kinematic fit to the 
hypothesis ■0(25') pir~pir + tt + ir~ tt + ir~ . The pir candidates for As are sent through a displaced vertex finder. Events 
which pass the fit with a fit probability greater than 0.01 are examined further. 

As the dominant decay mode in S* decays includes a A in the decay chain, a loose cut is placed on the pn" mass 
(|M p7r - — M\\ < 20 MeV/c 2 ) to enhance the signal fraction (Figure |l4|). Due to it combinatorics, each event that 
passes the kinematic fit is counted four times in this plot. 
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FIG. 14. Distribution of pn masses from a kinematic fit to ^(2S) 
S°(1530)H°(1530) candidate events. 



pp-K + -K 7r + 7r 7r + 7r showing the A peak in 



Similarly, as there is a S~ in the decay chain, a cut is made on the ptt~tt~ invariant mass; M pw - T - is required to 
be within 20 MeV/c 2 of the nominal mass of the S~, as shown in Figure 15. Due to ir combinatorics, each event is 
counted twice in this plot. 
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FIG. 15. Distribution of pn tv masses from a kinematic fit to ip{2S) — + ppTT + ir tt + -k tt + -k , showing the 3 peak in 
H(1530)°H(1530)° candidate events. 

All events which remain after the above cuts are graphed in Figure |l6| with Mp^+K+TT- on the vertical axis and 
M p7r - 7r - 7r + on the horizontal. The signal region is shown as a circle at (1.531,1.531). No events fall within the signal 
region defined as a 50 MeV/c 2 radius from the central value. We set an upper limit of 2.3 events at 90% CL for 



14 



2 



+ 



<M 

o 

> 

CD 

O 

C/D 
CO 
03 



CD 
CO 

■g 
c 

03 

O 

o 

* 

iW 



1.8 



1.6 




'1.4 1.6 1.8 2 

S*° Candidate Mass (GeV/c 2 ) 
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The circle denotes the signal region, 3 sigma from the nominal mass of the H(1530)°. 

d. 4>(2S) -* 

The dominant fl~ decay chain is Q~ — > AK~, A — > £>7r~ with a total branching fraction of 43% (?]]. We look for 
ip(2S) -> events with the topology -0(25") — » ppir + Tr K + K , i.e. six charged tracks where the j»7r and p7T + 

are consistent with being from the decay of a A or A. 

We select events with six charged tracks in the polar angle region |cos#| < 0.8 and with zero net charge. The 
remaining events are subjected to a 7-constraint kinematic fit to the hypothesis ip(2S) — > AAK + K~ , Max- = 
The fit is applied for each of the 36 particle assignment possibilities. Only the assignment with the best probability 
in the kinematic fit is considered. 

Figure |lj shows the AK~ mass distribution for the selected events, where the solid line histogram is data and the 
crosses are from Monte Carlo, normalized to three events. There are no candidates within three sigma of the nominal 
£!~ peak, thus an upper limit of 2.3 is assigned at the 90 percent confidence level. 
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V. CONTINUUM BACKGROUND 



A few percent of the hadronic events in our data sample originate from non-resonant e + e _ — > qq annihilation events. 
We use a 5.1 pb _1 data sample taken off resonance to determine the level of continuum e + e~ — » BiBi contamination 
to our event samples. We find no events that survive the analysis procedures and event selection criteria identical 
to those described above for either of the modes AA or A ++ A . We conclude that continuum events comprise a 
negligibly small contamination to our data samples. 



VI. DETERMINATION OF THE TOTAL NUMBER OF ip(2S) EVENTS 

We determine the number of tp(2S) events in our data samples from the observed number of cascade decays of 
the type tp(2S) — > tt + tt^ J/tp, J/tp — > X. The pions are reconstructed, and the recoil mass of the two pions is fit to 
determine the total number of tp(2S) — > tt + it~ J/tp events. From this fit, the total number of these events corrected 
for detection efficiency is 1.227 ± 0.003 ± 0.017 x 10 6 . The analysis for this is documented in reference jTTj]. 

The total number of tp(2S) events is determined by dividing the number of events in the previous paragraph by 
the PDG branching fraction for the mode ip(2S) — » tt + tt~ J/tp M. This number is determined to be 3.95 ± 0.36 x 10 6 
where the error is dominated by the error on the tp{2S) — > ir + n~ J/tp branching fraction EM. 



VII. ACCEPTANCE AND EFFICIENCY 



A. pp 



We determine the efficiency for tp(2S) — > pp events from a sample of Monte Carlo simulated events. Events were 
generated with a distribution of 

dN 2n 

oc 1 + a cos 



dcos 6 



with a = 0.61±0.23. This proportionality constant was measured in the J/tp — > pp system by the Mark II collaboration 
H and by the DM2 collaboration ||. Out of 20000 events generated, 14857 events survive these cuts, yielding a general 
efficiency of 0.743. The collinearity cut is also purely geometric, and has an efficiency of 0.999. 
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As the BES Monte Carlo is of limited usefulness for simulating detailed hadronic interactions, cuts which are affected 
by such must be corrected for by the examination of real data. Fortunately, there is a subset of events which allow 
the effects of th ese cuts to be determined. A clean sample of pp pairs was aquired from the analysis of A ++ A~" 
in Section IV A . The J/ip contamination shown in Figure || is the origin of this sample. These events are used to 
determine the Esc, and dE/dx cut efficiencies. 

This study is summarized in Table [| The systematic error was determined from both Monte Carlo statistics and 
variation of cuts. The product of all efficiencies is: 0.234 ± 0.022. 



Cut 


EMC 




Se 


523, eff 


SB, total 


General 


0.743 




0.006 


0.3 x 10" s 


0.3 x 10" s 


Muon ID 


0.696 




0.007 


0.2 x 10~ 5 


0.2 x 10" 5 


BSC Geom 


0.768 




0.009 


0.3 x 10~ 5 


0.3 x 10~ 5 


Esc < 0.7 




0.610 


0.048 


1.8 x 10~ s 


1.8 x 10~ s 


\XSP\ < 3 




0.968 


0.047 


1.1 x 10" 5 


1.2 x 10" 5 


Collinearity 


0.999 




0.001 


0.02 x 10~ 5 


0.02 x 10~ 5 



TABLE I. Relative efficiencies and systematic errors for 
cuts in the mode ip(2S) — > pp as modeled by Monte Carlo for 
geometric efficiencies and J/ip data for PID efficiencies. Last 
column includes combined systematic error due to variation 
of cuts. 



B. AA, E°E°, and E~E~ i 



We determine the efficiency for the hyperon-pair channels completely from Monte Carlo simulated events. Here we 
generated 20000 events in each mode with a 1 + a d cosd distribution^ = 0.67 ± 0.21, 0.22, 0.5 ± 0.5 for AA, 
and respectively. The value of aid was not varied for the S°E mode as the statistical error was large. The 

values for ay are those determined by the Mark II collaboration || and by the DM2 collaboration || . 

The resulting efficiencies are summarized in Table [□]. The systematic error reported is a combination of Monte 
Carlo statistics and variation of ay. Also, an additional 10 percent error is added because of uncertainties in the 
kinematic fitter used in these analyses. 

In these three modes, we require two A's that decay to charged pir final states, which have a branching fraction 
23(A — > p7r~) = 0.639 ±0.005" the other decay modes that are required are very nearly unity, namely B(E — ► A7) = 1.0 
and B(E~ — > A7r _ ) = 0.999 M. The branching fraction acceptance for each channel is the A — > pir~ branching fraction 
squared: 0.41 ±0.01. 



mode 


N evt 


efficiency 


B.F. Acceptance 


KFit 


pp 


201 ± 14 ± 21 


0.227 ±0.032 


1.00 




AA 


80 ± 9 ± 4 


0.27 ±0.01 


0.41 ±0.01 


10% 


E U E U 


8±3±2 


0.043 ± 0.003 


0.41 ±0.01 


10% 




12 ±3.4 ±0.2 


0.078 ±0.01 


0.41 ±0.01 


10% 


A++A— 


157 ± 13 ±34 


0.31 ±0.02 


1.0 ±0.01 


10% 


E*+E*" 


14±4±3 


0.104 ±0.005 


0.316 ±0.011 


10% 




< 2.3 


0.041 ±0.001 


0.172 ±0.001 


10% 




< 2.3 


0.042 ± 0.001 


0.187 ±0.001 


10% 



TABLE II. Number of events, efficiencies, branching fraction acceptances, additional systematic errors due to the kinematic 
fit for V(2S) BiBl. 
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C. Acceptance and Efficiency of the Decuplet Pairs 



We determine the efficiency for the decuplet hyperon-pair channels completely from Monte Carlo simulated events. 
Here we generated 20000 events in each mode with &l+ad cos 9 distribution, ad varying between and 1 for A ++ A~~, 
and constant at 0.6 for the other modes. The resulting efficiencies are summarized in Table ||, where the systematic 
error reported is a combination of Monte Carlo statistics and variation of ad- Also, an additional 10 percent error is 
added because of uncertainties in the kinematic fitter used in these analyses. 

The branching fraction for A ++ — > 7r + ir + is greater than 0.99 M, thus the branching fraction acceptance used is 
1.00 ± 0.01. The £*+£*- decay contains two As going to pir (0.63<F) and two E*s decaying to Att (0.88 2 ), for a total 
acceptance of 0.316±0.011. The H*°S*° decay has 3 components to the acceptance: S*° — »~ ir + (0.650 2 ), S~ — * Air~ 
(0.999 2 ) A — > pn~ (0.639 2 ), with a total acceptance of 0.172 ± 0.001, and r2~r2 + has only two components in the 
acceptance, Cl~ -> AK~ (0.678 2 ) and A -»■ pn^ (0.639 2 ), with a total acceptance of 0.187 ± 0.001. 



VIII. RESULTS 



The branching ratios B(ip(2S) — > BiBi)/B(ip(2S) — > J/tpir + Tr~) are calculated by dividing the number of events 
in each mode, corrected for efficiency and branching fraction acceptance, by the corrected number of events in the 
reference mode, as noted in Section [vj. The final branching fractions are determined by multiplying the above 
branching ratios by the PDG value for B(ijj(2S) — > J/i/i7r + 7r~), 0.310 ± 0.028. These are shown along with the 



branching ratios in Table III. 

We compare our results for the branching fractions to previous limits and results in Figure [l^. Our measured value 
for the B(ip(2S) — >£>p) is about one standard deviation higher than the previous DASP measurements, which was 
based on 4 events [|| and a Mark I measurement with similar statistics 0. The results for AA and are within 

the PDG upper limit values. There are no previous experimental results for E°E° or any of the decuplet modes. 



mode 


N evt , Corr 


B{B l B l )/B{J/ifjTi + -K-) (xHT 4 ) 


6(xl0~ 5 ) 


PP 


856 ±60 ±119 


6.98 ±.49 ±.97 


21.6 ± 1.5 ±3.6 


AA 


718 ± 80 ± 84 


5.85 ±.65 ±.69 


18.1 ±2.0 ±2.7 


E U E U 


456 ± 162 ± 152 


3.7 ± 1.3 ± 1.2 


12±4±4 




371 ± 108 ± 49 


3.0 ± .9± .4 


9.4 ± 2.7 ± 1.5 


A++A- 


506 ± 40 ± 127 


4.12 ±.33 ±1.04 


12.8 ± 1.0 ±3.4 


E*+E*~ 


419 ± 113 ±97 


3.4 ± .9± .8 


11±3±3 




< 322 


< 2.6 


< 8.1 




< 290 


< 2.4 


< 7.3 



TABLE III. Numbers of events corrected for efficiency and branching fraction acceptance, 
branching fraction B(ip(2S) -» B~B~) /B{ip{2S) -> J/^tt+tt") and final branching ratios for 
B(tp(2S) — > BiBi). Column 3 is calculated by dividing the corrected number of events in each 
mode by the corrected number of events in the reference mode. 
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FIG. 18. Comparison of measured branching fractions (circles) with previous measurements (triangles). Two previous 
measurements are upper limits. 

In Fig. [l9|, we plot the reduced branching fractions derived from our measurements. The results show a trend 
to smaller values for the higher massses, similar to that seen for the J/ip and are only marginally consistent with 
expectations from flavor-S'[/(3) symmetry. Higher precision measurements both for the J ftp and ip(2S) would clarify 
this issue. 




FIG. 19. The reduced branching fractions \Mi\ 2 = B(t[>(2S) -> BiBi)/(np* / yfs) for ip(2S) -> B t B, decays. 

A comparison to the perturbative QCD predictions of Bolz and Kroll Q is shown in Figure The results match 
quite well with these calculations. 
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FIG. 20. Comparison of B(ip(2S) — > BiBi) to Bolz and Kroll's predictions from perturbative QCD. Horizontal line is 

^exp/^qcd — 1.0. 



Our measured ip(2S) branching fractions agree with expectations derived from the application of the 12% rule to 
the PDG values for the corresponding J/tp decays for the modes pp, AA, S S° A ++ A and as 

shown in Table IV and in Figure ^ There are no results for J/ip — > S°(1530)S°(1530) and JJip — > is not 

kinematically allowed. 



Decay Mode 


B 


0.116 x B (xl0~ 5 )| 




(2.14 ±0.10) x 10 -3 


24.8 ± 1.2 


J/V>^AA 


(1.35 ±0.14) x 10 -3 


15.7 ± 1.7 




(1.3 ±0.2) x 10~ 3 


15. ± 2. 




(0.9 ±0.2) x 10~ 3 


10. ± 2. 




(1.10 ±0.29) x 10~ 3 


12.8 ±3. 




(1.03 ±0.13) x 10~ 3 


11.9 ±01.5 



TABLE IV. Branching ratio predictions for ip(2S). 
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FIG. 21. The ratio B(tp(2S) — » BiBi)/B(J/tp — » BiBi). Horizontal Line is the 12 percent ratio expected from factorizing 
the ip(2S) — > BiBi Feynman diagram. 



IX. CONCLUSIONS 



We report measurements of the branching fractions for ip(2S) — > pp, AA. 5 A ++ A and 

E + (1385)E~(1385), along with upper limits for the decays ip(2S) -> S°(1530)S°(1530) and fi - n+. The measured 
branching fractions agree with expectations based on an application of the 12% rule to the corresponding J/ip decays. 
The reduced branching fractions decrease with increasing baryon masses, showing some deviation from expectations 
based on Ravor-SU (3) symmetry. 
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